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Abstract 

We study the phase diagram and multicritical behavior of anisotropic Heisenberg antiferromag- 
nets on a square lattice in the presence of a magnetic field along the easy axis. We argue that, 
beside the Ising and XY critical lines, the phase diagram presents a first-order spin-flop line start- 
ing from T = 0, as in the three-dimensional case. By using field-theory methods, we show that 
the multicritical point where these transition lines meet cannot be 0(3) symmetric and occurs at 
finite temperature. We also predict how the critical temperature of the transition lines varies with 
the magnetic field and the uniaxial anisotropy in the limit of weak anisotropy. 

PACS numbers: 64.60.Kw, 05.10.Cc, 05.70.Jk, 75.10.Hk 
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I. INTRODUCTION 



Anisotropic antiferromagnets in an external magnetic field have been studied for a long 
time. In order to determine their phase diagram, they have often been modelled by using 
the classical XXZ model 

7~t = J } ] S n ■ S m + A ^ ^ S mtZ S njZ — H *y ^ S mjZ , (1) 

(mn) (mn) m 

where (mn) indicates a nearest-neighbor pair. Equivalently, one can use the Hamiltonian 

H = J Sn ■ S m + Dj2s 2 m , z -Hj^Sm^ (2) 

(mn) m m 

with a single-ion anisotropy term. The most interesting case corresponds to uniaxial systems 
that show a complex phase diagram. They correspond to Hamiltonians with A > or D < 0. 
Several quasi-two-dimensional uniaxial antiferromagnets have been studied experimentally, 
such as K 2 MnF 4 , Rb 2 MnF 4 , RbsMnCL;. 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 

Some general features of the phase diagram of anisotropic antiferromagnets are well 
known. Their phase diagram in the T-H plane presents two critical lines, belonging to the 
Ising and XY universality class, respectively, which meet at a multicritical point (MCP). 
The nature of the MCP has been the object of several theoretical studies. In three dimen- 
sions (3D), the issue has been recently studied using a field-theory approach. 8 ' 9 The starting 
point is the O(ni) © 0(n 2 ) symmetric Landau-Ginzburg- Wilson (LGW) $ 4 theory 10 

Hlgw = y"^{^[(^ 1 ) 2 + (^ 2 ) 2 ] + i[r 1 ^ + r 2 ^] 

+ ^Ko($?) 2 + «2o(^) 2 ] + ^o^}, (3) 

where $i and $ 2 are vector fields with n\ and n 2 real components, respectively. In our 
case nx — 2 and n 2 = 1 so that Hamiltonian (3) is symmetric under Z 2 © 0(2) transforma- 
tions. The RG flow has been studied by computing and analyzing high-order perturbative 
expansions, 8 ' 9 to five and six loops. It has been shown that the stable fixed point (FP) of 
the theory is the biconal FP, and that no enlargement of the symmetry to 0(3) must be 
asymptotically expected because the corresponding 0(3) FP is unstable, correcting earlier 
claims 10 based on low-order e-expansion calculations. The perturbative results allow us to 
predict that the transition at the MCP is either continuous and belongs to the biconal uni- 
versality class or that it is of first order — this occurs if the system is not in the attraction 
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FIG. 1: Phase diagram with a bicritical point where an Ising and an XY line meet. Here, T is the 
temperature and g a second relevant parameter. The thick black line ("flop line") corresponds to 
a first-order transition. 

domain of the stable biconal FP. Which of these two possibilities occurs is still an open issue. 
A mean-field analysis 10 ' 11 shows that the MCP is bicritical if 5q = Wio^2o — $ w o < 0; an d 
tetracritical if 5q > 0. Figs. 1 and 2 sketch the corresponding phase diagrams. At one loop 
in the e expansion the biconal FP is associated with a tetracritical phase diagram. 10 Since 
uniaxial antiferromagnets have a bicritical phase diagram, this result predicts a first-order 
MCP as in Fig. 3. First-order transitions are also expected along the critical lines that 
separate the disordered phase from the ordered ones. They start at the MCP, extend up to 
tricritical points, and are followed by lines of XY and Ising transitions. 

The nature of the MCP is even more controversial in two dimensions (2D) where different 
scenarios have been put forward; see, for example, Refs. 1,2,3,4,5,6,7,12,13,14,15,16. The 
existence of the Ising and the XY Kosterlitz-Thouless (KT) critical lines is well established. 
But, other features, like the existence of the spin- flop line, the nature of the MCP, and 
the MCP temperature, are still the object of debate. For example, the two different phase 
diagrams sketched in Fig. 4 have apparently been both supported by recent numerical Monte 
Carlo analyses: Ref. 15 claims that MC data are in agreement with the phase diagram 
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FIG. 2: Phase diagram with a tetracritical point where an Ising and an XY line meet. Here, T is 
the temperature and g a second relevant parameter. 

reported on the left, while Ref. 16 supports that reported on the right. Similar contradictions 
appear in the analysis of the experimental data. 1 ' 2,3 ' 4,5,6 ' 7 

In this paper we study the nature of the MCP in 2D within the field theoretical approach. 
Our main findings are the following: 

(i) The MCP cannot be 0(3) symmetric and cannot occur at T = 0. Indeed, the magnetic 
field and the anisotropy give rise to an infinite number of relevant perturbations of 
the 0(3) FP. Moreover, since at T = the order parameter is discontinuous as H 
is varied, we expect a first-order spin-flop line. We thus predict a finite-temperature 
MCP and exclude a phase diagram such as the one shown in Fig. 4 on the right. 

(ii) We study the limit of weak anisotropy. In this limit we determine how the critical 
temperature of the Ising and KT transition lines and of the MCP vary with H and A. 

Our results on the phase diagram of the classical XXZ model (1) are also relevant for 
the phase diagram of quantum spin-S" XXZ models 17 ' 18 and related systems, such as the 
hard-core boson Hubbard model. 19 ' 20 

We finally mention that anisotropic antiferromagnets modelled by the XXZ model (1) on 
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FIG. 3: Phase diagram with a bicritical point of first order. The thick black lines correspond to 
first-order transitions. 

a triangular lattice or on a stacked triangular lattice (and, more generally, on lattices that 
are not bipartite) are expected to show a different multicritical behavior, essentially because 
of frustration. Ref. 21 reports a study of the possible phase diagrams in the mean- field 
approximation and a field-theory study of the 3D renormalization-group (RG) flow. 

The paper is organized as follows. In Sec. II we summarize some general predictions for 
the isotropic model. In Sec. Ill we investigate the stability of the 0(3) FP by using field 
theory and give predictions for the critical temperature as a function of H and A for small 
A and H. Finally, in Sec. IV we discuss the implications of the field-theoretical results and 
discuss the possibile phase diagrams that are compatible with them. 

II. THE ISOTROPIC ANTIFERROMAGNET IN TWO DIMENSIONS 

For H = A = 0, the critical behavior of model (1) is well known. Indeed, if we perform 
the change of variables 
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FIG. 4: Phase diagrams in the T-H plane in the presence of easy-axis anisotropy A > in 
2D. The left figure shows a first-order transition line separating the antiferromagnetic (AF) and 
spin-flop (SF) phases, with a bicritical point where the Ising (Is) and the XY Kosterlitz-Thouless 
(KT) critical lines meet. In the right figure the Ising and KT lines meet at a zero-temperature 
0(3)-symmetric MCP. 

where, in 2D, [n] = n x +n y [n = (n x ,n y )], we obtain the ferromagnetic Heisenberg model 
for which several results are known. 22 In particular, if we define the two-point function 



G(n) = (S ■ S n ), 

a universal critical behavior is observed for the staggered susceptibility 

Xs = £(-l) lB] G(n) 

n 

and for the staggered second-moment correlation length 



c 2 = ^£(-i) w m 2 gw. 



(5) 



(6) 



(7) 



Under the mapping (4), Xs and £ s go over to the standard susceptibility and correlation 
length of the ferromagnetic model. In 2D, perturbation theory together with the RG allows 
one to derive the critical behavior of these two quantities for T — > 0. 23 Using the results of 



Refs. 24,25, we obtain 



Xs 

6 



C x i 2 s (J^) [1 + 0.18169T + 0.1334T 2 + 0.1346T 3 + 0(T 4 )] 
Qe 27r/T (L\ [i - 0.0914T - 0.1969T 2 + 0(T 3 )] . 



(8) 



The constants and C x cannot be determined in perturbation theory. Numerical values 
are reported in Ref. 26: 

C x = 93.25(3), Q = 0.0124783(12). (9) 

The asymptotic expansions (8) are quite accurate for T < 0.3, within a few percent at 
most. 27 ' 28 

III. CRITICAL AND MULTICRITICAL BEHAVIOR FOR SMALL A AND H 
A. General results 

We consider the 0(3) isotropic model at H — and add terms that break the 0(3) 
symmetry down to Z 2 <S> 0(2) (for instance the magnetic field or the anisotropy). The 
corresponding general Hamiltonian is 

n gcn = jJ2s n -s m + \Q(s). (io) 

(nm) 

If Q(S) is a relevant perturbation of the 0(3) FP, the 0(3) critical point is a MCP in the 
full theory. In 3D we can write the singular part of the free energy 10 for A — * as 

F sing ~u 2 t - a B(X), X = u x u; 4 >, (11) 

where a and are the 0(3) specific-heat and crossover exponents, respectively, B(X) is a 
universal scaling function, and u t and u\ are the scaling fields associated with the tempera- 
ture and with A. In general, we expect 

u t = t + k\ (12) 

where A; is a constant, t = T/T (z) — 1 is the reduced temperature, and T G ( 3 ) is the critical 
temperature of the isotropic model. No such mixing between t and A occurs in u\, since 
u\ vanishes for A = 0. Hence, we can take u\ — A. The crossover exponent is related 
to the RG dimension y\ of the operator Q that represents the perturbation of the MCP: 
(j) = y x v. Suppose now that the system has a critical transition for A ^ at T C (A). Since 
the singular part of the free energy close to a critical point behaves as (T — T c ) 2 ~ a , we must 
have B(X C ) = 0, where X c is the value of X obtained by setting T = T C (A). This equation 
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is solved by X c = X±, where X± are two constants that depend on the sign of A, such that 
X + > 0, X_ < 0. Hence, we obtain 



A 



T, 



0(3) 



= X ± . (13) 



It follows 



T C (A) = T 0(3) [l + (A/X + )^ -k\-\ ], A > 0, 

T C (A) = T 0(3 )[l + (X/X_)^ - fcA + • • •], A < 0. (14) 

These expressions provide the A dependence of the critical temperature for A small. Note 
that, depending on the sign of A, T C (A) varies differently. The sign of A may also be relevant 
for the nature of the phase transition. Indeed, the low-temperature phase may be different 
depending on this sign: in this case one observes critical behaviors belonging to different 
universality classes for A > and A < 0. 

One can generalize these considerations to the case in which there are two relevant per- 
turbations Q\ and Q 2 with parameters Ai and A 2 . In this case we can write 

^ sing ^u 2 t ' a B(X 1 ,X 2 ), (15) 

with 

X 1 = uiu^ 1 , X 2 = u 2 u^ <t>2 , 
u 1 — X 1 + c 1 X 2 , u 2 = Ai + c 2 A 2 , 

ut = t + fciAi + k 2 X 2) (16) 

where hi, k 2 , c±, and c 2 are constants. We have assumed here that there are two RG relevant 
operators that break the 0(3) invariance, with RG dimensions y\ > y 2 > 0, as is the case 
for anisotropic systems in the presence of a magnetic field. The crossover exponents are 
4>i — V V\ an d 02 = vyi- The scaling fields U\ and u 2 are linear combinations (and, beyond 
linear order, generic functions) of A x and A 2 , due to the fact that the lattice operators Q\ 
and Q 2 generically couple both RG operators. 

Let us now assume that, for Ai and X 2 small, the model shows generically two types of 
phase transitions belonging to different universality classes and that, for specific values of 
the ratio X 1 /X 2 , multicritical transitions occur. As before, we wish to derive the dependence 
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of the critical temperature as a function of Ai and A 2 for Ai, A 2 — > 0. Since 

X 2 = u 2 \u l \-^ /4 ' 1 \X 1 \^ / ' t, \ (17) 

we can rewrite Eq. (15) as 

•Fsing ~ u^B^X^u^^), (18) 

where we have introduced two different functions depending on the sign of u±. If u\ > the 
relevant function is B + (x,y), if u\ < one should consider B-(x,y). At the critical point 
we must have 

B±(X lc ,« 2 |« 1 |-^ 1 ) = 0, (19) 

which implies 

X lc = F ± (u 2 \ Ul \-^), (20) 



with F + (x) > and F-(x) < 0. We obtain finally 

i/4>i 

- fciAi - A; 2 A 2 . (21) 



T 

- 1 c 



0(3) 



Mi 



_F ± (u 2 |ui|~ < k / ^ 1 )_ 

Let us now discuss some limiting cases. First, assume that MalMil-^/* 1 < 1. This implies 
that X 2 is small compared to X 1: so that we can neglect X 2 . We are back to the case 
considered before and thus we can use Eq. (14). Depending on the sign of u x we have: 

1) If u\ > 0, the phase transition is located at 

'' 1 « ( ^rY* 1 ~ hXi - h\ 2: (22) 



^0(3) , 

where X + > 0. 
2) If ui < 0, there is a phase transition located at 



J o(3) 

where X_ < 0. 



' - 1 » f^-V M - *iAi - feAa, (23) 



Comparing with Eq. (21) we obtain F±(0) = _X±. 

The second interesting limiting case corresponds to ui — > 0. In this case X x can be 
neglected and we need to consider only X 2 . Therefore, depending on the sign of u 2 , we have 
a MCP located at 



1/4, 



2 



^- - 1 « f^- ) - fciAi - A; 2 A 2 , (24) 

'0(3) V A mc,±, 
9 



where X, 



mc,+ 



> and X 



mc. 



: j+ < 0. Consistency with Eq. (21) requires 



F±(x) « ±(x/X mC!+ )- 4 ' l/ ^ for x -> +oo, 
F±(ar) « ±(x/X mc _)-^^ 2 for x -> -oo. 



(25) 



One may devise simple interpolations that are exact for x — > and x — > oo. If, for instance, 
■ui > and -u 2 > 0, we may consider the approximate expression 



The functions F±(x) are crossover functions that interpolate between the two regimes in 
which only one of the relevant operators is present. As far as the nature of the transition, 
the relevant quantity is the sign of u±. In general, we expect that the transition belongs to 
different universality classes depending on the sign of u\. For u\ — the leading relevant 
operator decouples and thus we obtain a MCP whose nature may depend on the sign of U2- 
The previous results apply to the 3D model but cannot be used directly in 2D, since in 
this case T Q ^ = and v is not defined (£ s ~ e 27r / T for A = 0). To investigate the 2D case, let 
us consider again Hamiltonian (10), let us assume that Q(S) renormalizes multiplicatively 
under RG transformations and that its RG dimension is 2 with logarithmic corrections (as 
we shall see, this is the case of interest). The perturbative analysis of the scaling behavior 
of the free energy is analogous to that presented in Ref. 29. The correct scaling variable is 



where we used Eq. (8), and p is a power that can be computed by using the one- loop 
expression of the anomalous dimension of Q. Then, the singular part of the free energy is 
given by 




(26) 




(27) 



(28) 



The critical line is again characterized by Y = Y±, i.e. by 




(29) 



with Y + > and Y_ < 0. Solving this equation for T c , we obtain 



1 — p 




\ ln(V±/A) 
ln(y ± /A) 



+ ••• 



(30) 



ln(y ± /A) 
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The discussion is analogous in the case there are two relevant perturbations. We define two 
scaling variables Y\ and Y 2 as 



(31) 



and write the free energy as 

F sing = C 2 B(Y 1: Y 2 ). (32) 

If pi < p 2 , the critical behavior depends on the sign of uy. For u± > and u± < one 
obtains different critical behaviors. In the limit in which Y 2 is small and can be neglected, 
we can write 



47T 



\n(Y ± / Ul ) 



1-pi 



ln±ln(r ± / Ml ) 



(33) 



\n(Y ± / Ul ) 

where Y± are two constants such that Y + > and Y_ < 0. As observed in the 3D case, 
Eq. (33) [it is the analogue of Eqs. (22) and (23)] holds only if Y 2 can be neglected. Since 

P1-P2 



Eq. (33) holds only in the parameter region in which 



(34) 



u 2 < Mi ( ^\n(Y ± /u 1 ) 



P2-P1 



(35) 



i.e. far from the MCP at U\ = 0. For ^0 this region shrinks to zero and thus a correct 
formula requires the full crossover function. 
For u\ = 0, we have a MCP at 



™ 47T 

T c = 

ln(y mC)± /u 2 ) 



1 ~P2 



In \ ln(y mC)± /u 2 ) 
ln(F mCi± /M 2 ) 



(36) 



depending on the sign of u 2 . 



B. Effective Hamiltonians 



Let us now apply the results of the previous section to the XXZ model. For this purpose 
we consider the ferromagnetic Hamiltonian corresponding to (1) under the mapping (4), 



Hf = a n ■ a m - ^n,zO- m , z + H ^(-l) [n] <r„ iZ , 



(37) 



(mm) 
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where we have set J = 1 for convenience. We argue that Hamiltonian (37) is equivalent to 
the Hamiltonian 

H* = -J2°n-v m + J2 a ^ h > a )°^)i ( 38 ) 

(nm) 1>1 

where 2 i(cr) are the zero-momentum dimension- zero spin-2/ perturbations of the 0(3) FP, 
h = H/T, a = A/T (note that the partition function depends on h and a). The operators 
2 i(cr) can be constructed starting from the symmetric traceless operators of degree 21, see, 
e.g., Ref. 30. Explicitly, for I = 1,2 we have 



cM*) = £W,-i/3), 

n 

0,(a) = (< z - \o\ z + A) . 



(39) 



The functions a;(T, h, a) are smooth and vanish for h — > and a — > 0. It should be stressed 
that the two Hamiltonians are equivalent only for the computation of the leading critical 
behavior of long-distance quantities. 

The equivalence of TCf and H* for what concerns the critical behavior can be justified on 
the basis of the results of Ref. 10 obtained in the usual LGW approach. In the absence of 
anisotropy and magnetic field, i.e. for A = H = 0, model (1) is 0(3) invariant and thus its 
critical behavior is described by the usual LGW Hamiltonian 



(40) 



where $ is a three-component vector. The magnetic field and the anisotropy break the 0(3) 
invariance. If we define $ = ((f), (p), where is a two-component vector and (p a scalar, the 
LGW Hamiltonian (3) corresponding to (1) can be written as 10 



H = H&+ / d d x 



where 



/22-P22 + /42-P42 + /44-P44) 



1 



(41) 



22 



p 2 - -$ 2 , 

^ 3 ' 



-P42 — $ -P22, 

r> ,o 4 ^<f> 2 ,o 2 1 ^ ^2^2 

F 44 = Lp - -<P if + — (<P ) , 



(42) 
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and fij are coupling constants depending on a, h 2 , and T. We have introduced here the 
homogeneous polynomial P m \. The polynomial Pu has degree / and transforms irreducibly 
as a spin-/ representation of the 0(3) group. 31 Polynomials P m i, m > I, are defined as 
P m i = ($ 2 )( m ~0/2p^ The classification of the zero-momentum perturbations in terms of 
spin values is particularly convenient, since polynomials with different spin do not mix under 
RG transformations and the RG dimensions y m i do not depend on the particular component 
of the spin-/ representation. We refer to Ref. 8 for details. Thus, in the LGW approach the 
0(3) Hamiltonian is perturbed by spin-2 and spin-4 perturbations. If one considers higher 
power of the fields, also spin-6, spin-8, . . . perturbations appear. They are irrelevant in 
3D, but must be considered in 2D, since in this case the field <3> is dimensionless and any 
polynomial of the fields is relevant. 

C. Three-dimensional behavior 

In 3D the critical behavior is well known: 

(i) In the absence of anisotropy, i.e. for A = 0, there is an XY critical line Txy(H) ending 
at the 0(3) critical point for H = 0. The phase diagram is symmetric under H — > —H . 

(ii) In the absence of magnetic field, i.e. for H = 0, there are two critical lines: an XY 
critical line for A < and an Ising critical line for A > 0, which meet at the 0(3) 
critical point as A — > 0. 

(iii) If A and H are both present, one observes two different phase diagrams depending on 
the sign of A. If A is negative, there is an XY transition. If A is positive, there is an 
Ising critical line for H small and an XY critical line for H large; the two lines meet 
at a MCP corresponding to a nonvanishing value of H. 

We wish now to determine the dependence of the critical temperatures on A and H close 
to the 0(3) point. For this purpose we investigate the relevance of the perturbations. In 
3D we expect that only polynomials of degree two or four are relevant. Thus, P22, P42 
and P44 are the only quantities that should be considered. The RG analysis is reported in 
Ref. 8 and indicates that P42 is irrelevant, while P22 and P44 are relevant. Perturbative field- 
theory calculations provide estimates of the corresponding RG dimensions: 8,32 y 2 2 ~ 1-79, 
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FIG. 5: Phase diagram of isotropic antiferromagnets in the T-H plane. 

?/44 ~ 0.01, and ~ —0.55. The relevant scaling fields are u\ = a + c\h 2 associated with 
P22 and u 2 = a + c 2 h 2 associated with P 44 . The scaling fields depend on h 2 , because of the 
symmetry under H — > — H transformations. Since P22 is the most relevant operator, the 
critical behavior depends on the sign of u±. We obtain XY behavior for a + c\h 2 < and 
Ising behavior for a + c\h 2 > 0. Since for A = only the XY transition is observed, the 
constant c\ must be negative. Hence, Ising behavior can only be observed for a positive and 
h 2 not too large, in agreement with experiments. The critical temperatures behave as 

T c ,xy - To(3) = b XY (-a - Cl h 2 ) 1 '^ - ha - k 2 h 2 , 
T c ,i s - T 0(3) = b ls (a + Cl h 2 yl^ - ha - k 2 h 2 , (43) 

where bxy, h s , h and h are constants, with bxv > 0, &i s > 0. Here 8 2 2 = 1.260(11), 
1/022 = 0.794(7). 

Multicritical behavior is observed for a = —c\h 2 (since c\ is negative this equality can 
only occur for a > 0). Since 1/044 ^ 100, the leading nonanalytic term with exponent 1/^44 
cannot be observed in practice, and therefore 

P c , mc - T (3) = -ha - k 2 h 2 . (44) 
D. Two-dimensional behavior 

The critical behavior in 2D is analogous to that observed in 3D. For A < there is an XY 
Kosterlitz-Thouless (KT) transition for any H, while for A > there is an Ising transition 
for small \H\, an XY transition for large \H\, and a MCP in between. The phase diagram 
in the two limiting cases A = and H = is reported in Figs. 5 and 6. 
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FIG. 6: Phase diagram of anisotropic antiferromagnets in the T-A plane for H = 0. 

To compute the position of the critical point, we should investigate the relevance of the 
perturbations P m [. As shown in Ref. 33, in 2D any spin-/ perturbation Pu is relevant at the 
0(3) FP, since the corresponding RG dimension yu is positive, indeed yu = 2 apart from 
logarithms which can be computed by using perturbation theory. One can also argue that 
perturbations P m i with m > I can be neglected. Indeed, spin waves, that are rotations of 
the spins, are the critical modes of the ferromagnetic 0(3) model. Changes in the size of 
the field $ should not be critical. Therefore, P m i = (cl? 2 )*™ - ')/ 2 ^; should be equivalent to 

($2)(m-l)/2 Pij _ where $ 2 ig the 

average of $ 2 . Thus, one should only consider the 
operators Pu. Equation (38) then follows immediately. The equivalence of model (1) with 
Hamiltonian (38) was already conjectured in Ref. 34, even though there only the leading 
spin-2 term was explicitly considered. 
The relevant scaling variables are 

* = «**"■. (46) 

which are associated with the spin-2/ perturbation. The power of T, which is universal, has 
been determined by using the perturbative results of Ref. 33, which provide the anomalous 
dimension of any zero-dimension spin-/ perturbation of the 2D iV-vector model. The scaling 
fields ui are linear combinations of the parameters that break the 0(3) invariance. As before, 
we write them as u\ = a + q/i 2 . As a consequence, the free energy can be written as 

^sin g ~r 2 ^,n,...)- (46) 

Note that different powers of T appear in the definition (45). The most relevant term for 
T — > corresponds to the spin-2 operator, since Y 2 i ~ Y 2 T l ^ 2l+1 ^ 3 and Z(2/ + 1) — 3 is 
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positive for I > 2. Next one should consider the spin-4 scaling variable Y 4 . If one considers 
the scaling limit at fixed Y 2 ^ or Y± ^ 0, the higher-order spin variables go to zero as 
T — > and thus represent corrections to scaling proportional to powers of T. Thus, in the 
scaling limit one can neglect Y e , Y 8 , . . ., and write 



^Sing~r 2 ^2(l2,l4), 

which coincides with Eq. (32). We can thus use the above-obtained results. 

(i) For ui < 0, i.e. a + cih 2 < 0, there is a KT transition at 

.1h±1h(Yxy/«i) 



(47) 



T, 



c,XY 



Aix 

ln(y XY /«i) 



1-5- 



(48) 



ln(y XY /«i) 

where Y X y < 0. Note that the occurence of an XY transition for a = implies c\ < 
as in the 3D case. 



(ii) For ui > 0, i.e. a + Citi 2 > 0, there is an Ising transition at 

47T 



T c ,\s — 



\n(Y ls /ui 



, \n\\n{Y ls / Ul ) 



(49) 



where F Is > 0. Since C\ < this transition can only occur for a > and h 2 small, 
(iii) For u 1 = 0, i.e. a + Cih 2 = 0, there is a MCP at 

47T 



ln(F mc //i 2 ) 



i _ i2 \n\HY mc /h 2 ) 



(50) 



ln(y mc //i2) 

where F mc > 0. Note that since the MCP occurs for a = —Cih 2 , we can replace 
U2 = a + c 2 /i 2 simply with h 2 . 

Note that, as discussed before, Eqs. (48) and (49) are valid only as long as condition (35) is 
satisfied, i.e. far from the MCP u\ = 0. 



IV. DISCUSSION AND CONCLUSIONS 

In Sec. HID we have discussed the behavior of the classical 2D XXZ model (1) on a 
square lattice in the presence of a magnetic field along the easy axis, for A and H small, 
focusing on the dependence of the critical temperature on A and H. In this section we wish 
to discuss the possible scenarios for the nature of the MCP. 
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T 

FIG. 7: Phase diagram in the T-H plane in the presence of uniaxial anisotropy A > 0. 

First, we exclude the possibility that the MCP is 0(3) symmetric and located at T = 0, 
and thus we exclude a phase diagram such as the one shown on the right in Fig. 4. This 
is implicit in the results of Sec. Ill D, but we wish to present here an extended discussion. 
There are essentially two observations that exclude an 0(3) MCP. 

• For T = the XXZ model (1) on a square lattice can be easily solved. 12 ' 35 At fixed 
A > 0, one finds two critical values of the magnetic field: 

H cl = AJV2A + A 2 , 

H c2 = 4J(2 + A). (51) 

For \H\ < H c i the system is in a fully aligned antiferromagnetic configuration, for 
H c i < \H\ < H c2 the system is in a spin-flop configuration, while for \H\ > H c2 all 
spins are aligned with the magnetic field. The spin- flop transition at \H\ = H c i is of 
first order, since the order parameter, the staggered magnetization, has a discontinuity. 
If the transition were at T = 0, the MCP should coincide with the spin- flop point 
| if | = H&. Thus, the magnetization, the susceptibility and all critical quantities 
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FIG. 8: Some possible phase diagrams inside the blob appearing in Fig. 7. 

would have a discontinuity at T = when varying H. It is unclear how an 0(3) 
critical behavior might be consistent with this discontinuity. On the other hand, the 
first-order behavior at T = is consistent with the existence of a first-order spin-flop 
line. 

• The second objection is based on the LGW analysis presented in Sec. Ill B. The 0(3) 
FP is unstable under an infinite number of perturbations. Thus, an infinite number 
of tunings is needed to recover the 0(3) symmetry. 

On the basis of these remarks, in 2D anisotropic antiferromagnets in a magnetic field along 
the easy axis, a scenario based on a T = 0(3)-symmetric MCP appears untenable. These 
conclusions are analogous to those that hold in 3D. 8,9 It should be noted that in 2D the 
argument is much stronger. While in 3D these conclusions are based on the numerical 
determination of the RG dimension of the spin-4 perturbation, in 2D they follow from the 
relevance of the spin-/ perturbations, which is an exact result. Moreover, while in 3D the 
0(3) behavior can be obtained by tuning a single additional parameter in such a way to 
decouple both the spin-2 and the spin-4 operator, in 2D the 0(3) behavior can never be 
obtained at finite A, since an infinite number of tunings is needed. 

Since we have proved that in 2D the MCP cannot have 0(3) symmetry, an open question 
concerns the nature of the MCP. We shall now show that the decoupled FP, corresponding 
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to a multicritical behavior in which the two order parameters are effectively uncoupled, is 
stable. The stability of the decoupled FP can be proved by nonperturbative arguments. 36 
Indeed, the RG dimension y w of the operator <J>^$| that couples the two order parameters 
in the LGW theory (3) is given by 

y w = — + — -2 = -l<0. (52) 

^KT Vis 

since z/kt = oo and u\ s = 1. Therefore, the perturbation is irrelevant and the decoupled FP 
is stable. Note that a decoupled MCP is always tetracritical, as in Fig. 2. 

The stability of the decoupled FP and the instability of the 0(3) FP is also consistent 
with some general arguments. 37 At a MCP the exponent 77 describing the critical behavior of 
the correlation function of the order parameter is replaced by a matrix 77^, see, e.g., Ref. 10. 
The conjecture of Ref. 37 states that Trr/ should have a maximum at the stable FP. This 
indeed occurs in the present case: at the decoupled FP we have Trr/ = 3/4 while at the 
0(3) FP Trr/ = 0. 

In Fig. 7 we show a plausible phase diagram, with three transition lines: a spin-flop first- 
order transition line, an Ising and a KT critical line. The phase diagram inside the blob is 
an open issue. Some possibilities are shown in Fig. 8: 

• Fig. 8 (a) presents a bicritical point, which may be associated with a stable biconal 
FP, whose attraction domain is in the bicritical region of the bare parameters of the 
LGW $ 4 theory. We should say that we do not have any evidence for the existence of 
such a FP. 

• In Fig. 8 (b) we show a tetracritical point, where the Ising and KT lines intersect each 
other. In this case the MCP may be controlled by the stable decoupled FP discussed 
above. 

• In Fig. 8 (c) the transition at the MCP is of first order. Starting from the MCP, 
the first-order transitions extend up to tricritical points, where the Ising and the KT 
critical lines start. This occurs if the system is outside the attraction domain of the 
stable FP of the RG flow. This scenario resembles the one predicted in 3D, see Fig. 3. 

• Finally, we cannot exclude phase diagrams like those shown in Fig. 8 (d) and Fig. 8 
(e). Case (d) is apparently observed in antiferromagnets with single-ion anisotropy 

19 



and more than nearest-neighbor interactions, 38 and also in hard-core boson systems, 20 
which are equivalent to anisotropic spin-1/2 XXZ systems in a magnetic field. 

Of course, further experimental and theoretical investigations are called for to conclusively 
settle this issue. 
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